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The composite operator formalism is applied to QCD at finite temperature 



to calculate the masses of scalar and pseudoscalar mesons. In particular the 
ratio of the sigma mass to the pion mass is an interesting measure of the 
degree of chiral symmetry breaking at different temperatures. We calculate 

X 

the temperature T* at which M a (T) < 2M 7r (T), above which the sigma partial 
width into two pions vanishes. We find T* = 0.95 T c (where T c is the critical 
temperature for the chiral phase transition), within the full effective potential 
given by the formalism. We find that an expansion a-la Landau of the effective 
potential around the critical point in the limit of small quark mass provides 
for a very good determination of T* . 
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I. INTRODUCTION 



The concept of chiral symmetry breaking in strong interactions goes back to the pio- 
neering work by Nambu and it has been widely explored and discussed since that time. 
Restoration of chiral symmetry at high densities was proposed by Lee and Wick |2|] already in 
1974. Theoretical procedures for studying chiral symmetry restoration at high temperatures 
were proposed on the same year by Dolan and Jackiw || and by Weinberg M. 

QCD at finite temperature and density has attracted much interest. For latest studies 
of the QCD phase diagram see for instance 

Physical applications are to high energy heavy ion collisions || and to the physics of the 
early universe 0. 

A review of lattice calculations can be found for instance in Ref . || . The Nambu- Jona- 
Lasinio model has been used under different approximations to obtain indications on the 
chiral transition (see for instance [p|-|14|j ) . 

A number of peculiar features may emerge due to the chiral phase transition and re- 
cently much attention has been attracted by the possibility of disoriented chiral condensate 
formation (see for instance |]l7)|-|T7|j ) . 

It is strongly suspected that a single transition occurs rather than separate transitions 
for deconfinement and chiral-symmetry. We had suggested a heuristic argument indicating 
that, at least for zero density, the critical temperature for chiral transition, T c , coincides with 
that for deconfinement, Td ||18|| . The order parameters usually considered cover extreme and 
opposite ranges. The thermally averaged Polyakov loop is suitable in the limit of infinite 
quark masses to describe the transition from the confined to the non-confined phase. The 
other extreme is the limit of vanishing quark masses, where the thermally averaged quark- 
antiquark bilinears are the typical order parameters for chiral symmetry transition. 

We shall deal here with the chiral transition, concentrating on the thermally averaged 
quark bilinears at finite temperatures as order parameters for chirality. For light massive 
quarks, such as u and d, rigorously, chiral symmetry is already broken in the Lagrangian, 
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but we can still retain the notion of phase transition, by looking at the region of T where the 
condensate has a rapid variation. The current quark mass plays a role analogous to that of 
an external magnetic field in the ferromagnetic transition, as it explicitly violates the chiral 
symmetry whose restoration characterizes the phase transition. 

The analysis will be based on a composite operator formalism at finite temperature. 
The formalism makes use of an effective action for composite operators [^,^]. Actually 
the effective potential admits of a Landau expansion around the critical point and thus the 
behaviour of the condensate is well reconstructed by knowing the coefficients (which are 
infrared safe) and the critical exponents. We shall compare with the results following from 
a Landau expansion and show the general agreement. 

Within the composite operator formalism we had already discussed the T dependence 
of f n {T) in the whole range of temperatures up to T c . We shall here mainly deal with 
the ratio of the scalar to pseudoscalar mass M^/M^ at varying temperature. This ratio 
has a peculiar theoretical interest as a sensible indicator of the degree of chiral symmetry 
breaking. 

Let us for the moment neglect the small u and d quark masses. When chiral symmetry 
is restored the two mesons are degenerate in mass. The ratio is then equal to unity. At zero 
temperature instead in the broken chiral phase the pion is a goldstone and has vanishing 
mass, while the sigma has a finite mass from the chiral condensation. Thus one expects the 
mass ratio to decrease from oo to one at the chiral transition. Quark masses will change 
this picture quantitatively leading to a decrease from a finite value to a value close to one. 

Beyond some temperature, before approaching the critical value for the transition, the 
sigma will not have phase space left to decay into two pions. Its instability will decrease 
when increasing the temperature, and the channel into two pions will finally be suppressed 
at some temperature in the vicinity but lower than the critical temperature. The decay 
channel into two photons will still be available and become the dominant channel |T3| , p2| . 

It is too early to say whether a possible experimental signature for the transition may 
be inferred from such behaviour. An accurate determination of the mass ratio versus tem- 

3 



perature is in itself of theoretical interest as it constitutes a significant parameter for the 
degree of chiral symmetry breaking at finite temperature. This we shall do in this note in 
the composite operator formalism by making use of parameter values obtained from fits to 
zero temperature QCD. 

In Section 2 we shall summarize the main results of the application of the composite 
operator formalism in QCD. In Section 3 we determine the relevant QCD parameters and 
discuss the hadronic masses at non zero temperature. The temperature dependence of the 
observables is discussed in Section 4, where we also calculate the temperature at which the 
sigma can no longer decay via strong interaction. 



II. EFFECTIVE ACTION FOR QCD COMPOSITE OPERATORS 



Following Ref. [0 the zero temperature Euclidean effective action for an SU(N) QCD- 
like gauge theory is 

5T 2 



rs 



-Tr In 



S" 1 



5S 



— Tr [SS] — T2(S) + counterterms 



(1) 



where Sq 1 = (ip — m), m is the bare quark mass matrix and r 2 (S) is the sum of all the two- 
particle irreducible vacuum diagrams with fermionic propagator S and S = — 5T 2 /SS. At 
two-loop level T 2 = |Tr(SAS), where A is the gauge boson propagator, so that £ = —AS, 
Tr [S 5r 2 /5S] = 2T 2 and one can rewrite Eq. ([!]) in terms of S 



T(S) = -Trln 
= — Tr In 



S 1 - S 



S" 1 



+ T 2 (S) + counterterms 



-Tr (SA^s) + counterterms 



(2) 



Here the variable 5] plays the role of a dynamical variable. At the minimum of the 
functional action, that is when the Schwinger-Dyson equation is satisfied, S is nothing but 
the fermion self-energy. A parametrization for I], employed in [TJJ was 



(s + «7 5 p)/(fc) = S s + r/5£ f 



(3) 
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with a suitable Ansatz for f(k), and with s and p scalar and pseudoscalar constant fields 
respectively to be taken as the variational parameters. 

The effective potential one obtains from Eq. (|2|) (see Ref. ]19|) is 

t, r 8 ^ 2 N r d 4 k r ^ _ _ 

V= — = = / - — -rtr £ s D fc £ s + S D D fc SJ - 

Q, 3C 2 g 2 J (2tt) 4 L s k s v k v \ 

-N Tr In \ik - (m + S s ) - + 5Z tr(m s) (4) 

where C 2 is the quadratic Casimir of the fermion representation (for SU(3) C C 2 = 4/3). 
Furthermore E s = X a s a f(k) /y/2, S p = X a p a f(k)/ y/2, m = X a m a /y/2 (a = 0, • • ■ , 8, 



A = y2/3 /, Aj = Gell-Mann matrices, % — 1, • • • , 8), g is the gauge coupling constant 
and Q is the four-dimensional volume. In Eq. (^) has a divergent part to compensate 
the leading divergence proportional to Tr(m s) in the logarithmic term. We remark that 
both the current mass m and the self-energy S are, in general, matrices in flavour space. 



However, as discussed in Refs. |]19| , |20| if we neglect the mixing between different flavours 



originating, for instance, from terms such as the 't Hooft determinant p3|, it follows that 
only the flavour diagonal elements of the fermion self-energy and mass can be different from 
zero at the minimum. With vanishing off-diagonal terms, the effective potential decomposes 
into the sum of n/ contributions (71/ = number of flavours), one for each flavour. Therefore, 
to study the minima of the effective potential, it is formally sufficient to consider a single 
contribution. Of course, the choice of a given flavour number will reflect in the particular 
parameters assumed. In the present paper, as in Ref. |19[], we will take rif = 3 and a number 
of colour N = 3. The value of the parameters will be specified later on. 

As far as the choice for the function f(k) is concerned, in QCD, the operator prod- 
uct expansion suggests (neglecting logarithmic corrections) to take for f(k) a momentum 
behaviour as 1/k 2 for large k 2 . We have chosen as a variational Ansatz |24j 



/(*) = pfss (5) 

where M is a momentum scale which is expected to be of the order of Aqcd- 

To extend the zero-temperature theory to finite temperature we can still work with 
continuous energies by substituting for the sum over discrete energies u n = (2n + 1)tt//3 



(where (3 = 1/T) a sum of integrals over continuous energies by means of the Poisson's 
formula |25| 

d A k r rfSfc 
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This substitution corresponds to the imaginary-time formalism after using the Poisson's 
formula and allows for an equivalent version af the Dolan- Jackiw finite temperature Feynman 
rules ii. 

In conclusion the final form for the effective potential, for a quark of mass m, is (see Ref. 
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where E s = s /(A;) and E p = p /(&)• 

Let us now comment on the choice for the gauge coupling constant. As suggested by 
asymptotic freedom and renormalization-group considerations, we expect the strong forces to 
weaken at high temperature |26|| . We shall then assume that in the UV region the coupling 
constant g(T) depends logarithmically on the temperature T. We take into account this 
assumption by writing 

g\T) _ l _ l 



2tt 2 c(T) Cq + c^T) 



(8) 



where b = 24it 2 /(11N — 2n/) and we will discuss the parameters Co, M and £ later on. 

At T ^ the effective potential does not acquire any extra divergence with respect to 
the T = case. The renormalization at T = can be performed by adding a counterterm 
and requiring that the derivative of the effective potential with respect to the term which 



breaks explicitly the chiral symmetry, evaluated at the minimum, satisfies for each flavour 
the renormalization condition [[HJ] 

5V 



lim ■ 



(9) 



The fermion condensate (V'V')t is related to the minimum s(T) of the effective potential, 
renormalized at the scale M through the relation (see Ref. 



W3 

(Wt = s(T) (10) 

and (^)o = (V»V)t=o- 

At finite T we do not have any additional divergence with respect to the T = case; 
nevertheless, in order to satisfy the generalization of the normalization condition (|5|) at finite 
T we have to add a finite counterterm to that determined at T = [27||. Finally we recall 



that with the appropriate normalization for the pion field (see next Section), Eq. (§) and its 
generalization at finite T is also equivalent to the Adler-Dashen formula 

Ml{T)fl{T) = -2m(#) T . (11) 



III. HADRONIC MASSES ATT/0 AND DETERMINATION OF THE 

PARAMETERS 

To compute the masses of the scalar and pseudoscalar mesons (a and n respectively) 
one has to take the second derivative of the effective potential (|7[) with respect to the scalar 
field s and the pseudoscalar field p, evaluated at the minimum (for the sake of simplicity we 
are assuming that the up and down quarks are degenerate in mass and m = {m u + rrid)/2 
in Eq. ©). 

The actual values of the masses will be obtained by multiplying the second derivative 
by the appropriate factor a that relates the physical fields ip n ((p a ) to p (s) according to 
ip n = a p, ip a = as. This factor can be obtained in terms of the decay constant f n 



through standard arguments of current algebra [pLOfl . One gets 
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(12) 



where s is the extremum of the effective potential in the presence of a bare mass. 

To derive a more physically transparent expression for the masses of the mesons as 
expressed by Eq. ([12]) through the second derivative of the effective potential, we will use 
the gap equation and the generalization at finite T of the normalization condition (|9|). The 
extremum condition is 
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+ m 5Z = (13) 
(14) 



where S a (fc) = s f(k). 

As far as the normalization condition is concerned, by using the relation (|T0|) between 
the scalar field at the minimum and the scalar condensate we shall write in general 

d 4 k S (fc) JnPko 



5Z = N 



M 3 4 

S n 



(15) 



2^' ' s ^ ' J (2ir)*k 2 + X 2 (k) 
with Eq(A;) = Sof(k) and where Sq is the minimum of the effective potential in the massless 



case. 



Now, by using the gap equations (p^3])-(|I4T) and the normalization condition ( |T5| ) in the 
explicit expressions ( |i~2|) for the meson masses, we can eliminate 5Z and c(T) to finally 
obtain, for a quark of mass m: 
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We notice that Eq. (|T6|) reproduces, in the small mass limit, the Adler-Dashen formula 
(PI) and differs from this formula only by terms of order m 2 . 

In order to fix the values of the parameters Co and sq(T = 0) we follow the same pro- 



cedure of Ref. ||19|| . In the massive case the effective potential is UV divergent and from 
the normalization condition (in the small-mass limit) and gap equation one is able to fix, 
through a self-consistency relation, the values Cq = 0.554. By inserting this value in the gap 
equation for massless quarks, one finds that chiral symmetry is spontaneously broken (at 
T = 0) via a minimum of V located at sq(T = 0) = —4.06 and the point sq(T = 0) = 
is a local maximum. To determine the mass scale M and the quark masses from the ex- 
perimental data, one has to derive the explicit expressions for the masses and for the decay 
constants of the pseudoscalar octet mesons which are the pseudo-Goldstone bosons of chiral 
symmetry breaking. These expressions constitute a system of coupled equations which we 
have solved by an approximation method. The experimental inputs are the decay constant 
and mass of the charged pion, the charged kaon mass and the electromagnetic mass differ- 
ence between the neutral and the charged kaon. The outputs of the numerical fit for the 
octet meson masses (agreement within 3%) are the masses of the u, d and s quarks and the 



mass parameter M. The values we get are the following [|19p0"| , p8 

M = 280 MeV 
m u = 8 MeV 

m d = 11 MeV (18) 

m s = 181 MeV 

m u + m d 

m = = 9.5 MeV 

2 

Finally, the parameter £, which appears in the expression (|S|) for the running coupling 
constant, has to be determined on phenomenological grounds. If we compare our model in 



the low-T regime with the results of Ref. |2{| we find £ = 0A4(rij — l)/n/ for rif flavours, 
which gives £ — 1 for rif = 3 

With these values for the relevant parameters it is easy to determine the masses of the 
pseudoscalar and scalar mesons at T = (we are using the value f n = 93 MeV for the pion 
decay constant) 

M n ~ 135 MeV, M a ~ 630 Mel/ (19) 

and (again at T = 0) the fermionic condensates at M for massless and massive quarks 

(^V')o — (—197 MeV") 3 (massless quarks) (20) 
(^V')o — (—200 MeV) 3 (massive quarks). (21) 

The value we have obtained for the mass of the scalar meson, M a ~ 630 MeV, accords 
very well with those obtained by other authors within different approaches to the study of 



dynamical chiral symmetry breaking in QCD 113,14 . 



IV. THE SCALAR-PSEUDOSCALAR M a (T)/M n (T) MASS RATIO 

It is well known that as the temperature T increases we expect some changes in the 
hadronic properties of matter as a consequence of a possible phase transition. 



For instance, a physical quantity as f^(T) will decrease as T increases pl |. This means 
that the pion tends to decouple from matter (quarks and leptons) as T increases. As far 
as the pion mass is concerned, M V (T) has a weak dependence on the temperature. In fact 
the value of M 7r (T) is dominated by the current quark mass for temperatures below the 
critical value, whereas it becomes independent of m at the critical temperature or above it. 
This should be expected because, for T ~ T c , the pion becomes an ordinary resonance and 
its mass is dominated by Aqqd and not by m. This is a clear sign that the pion loses its 
Goldstone nature once the approximate chiral symmetry is restored ||20| . 

The scalar meson mass M a (T) will decrease in association with partial chiral symmetry 
restoration in hot matter. In fact the scalar a meson can be described as a qq quasi bound 
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state and its mass is about twice the constituent quark mass [ 13| , pL9 , |30 1 . Thus the a-mass will 
decrease substantially with increasing temperature and this behaviour is mainly determined 
by the T-dependence of the constituent quark mass. The decreasing of the a-mass is expected 
from the restoration of chiral symmetry at T = T c . In fact M a — > for m = at T = T c . 
As a consequence M a (T) will become smaller than 2M n (T) at some temperature T* < T c 
and the cx-meson, which has a large width at T = due to the decay o — > 2tt, would appear 
as a sharp resonance at high temperature since the amplitude for a — > 2ir would vanish 
above T* (near the critical point of the chiral phase transition). Thus one can expect to 
have a better chance of seeing the a-meson emerge as a clear, narrow state at non zero 
temperature. In fact, since for T > T* the decay width of o by strong interaction vanishes, 
only the electromagnetic width coming from a — > 27 is left ||L3| ) p2|| . 

Within our model we are able to make some rough estimate of the temperature T* at 
which this phenomenon takes place by simply studying the behaviour of the meson masses 
(|T6|), (|I~7|) as a function of the temperature. However, since we are, as already said, mainly 
interested in determining the temperature T* at which M cr (T) < 2M 7r (T), instead of evalu- 
ating separately M a {T) and M n (T) it will be sufficient to study the ratio M^/M| 

Ml ( 24 

Ml ~ m 
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By straightforward calculations at the leading order in a = m/M we obtain 
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(23) 



Furthermore, as we expect that the temperature T* will be not far from the critical tempera- 
ture of the chiral phase transition, we can start our investigation by performing a mean-field 
expansion a-la Landau of the effective potential in the small quark-mass limit around the 



critical point. Following Ref. 20] we have 
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V 
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a 2 (T) <^ 2 + a 4 (T) + 61 (T) a s + 



(24) 



where (p 2 = s 2 + p 2 . 

If we now compute the masses of the scalar and pseudoscalar mesons from the expression 
of the effective potential we get 



M 2 



1 



' M ° 4(T) Sq(T) + ... 



(25) 



M2 m h(T) 

By noticing that the generalization at finite T of the normalization condition @ on the 
effective potential reads now 



b x (T) = 2 c(T) 

and that the coefficient a 4 (T) is given by the expression (see Ref. JH] 

A 1 



(26) 



a 4 (T) = 4vr 2 



inj3My 
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J/ 2 (l/ 2 + I) 2 + s 2 



(27) 



it is easy to check that Eq. (|25|) is nothing but Eq. (|23|) . 

Just to have a rough idea of the value of the temperature T* we can push our approxi- 
mation further by evaluating a 4 (T) and b\ (T) at T c and using for So(T) the expression for 
the massless minimum around T c p0| 

j 1 \ 1/2 



s (T) = a s (T c ) 1 



T < T, 



(28) 



in Eq. 

The numerical values for the coefficients relevant for the evaluation of the mass ratio are 



a 4 (T c ) = 0.0033 , h(T c ) = 1.25 , a s {T c ) = -7.28 



(29) 



and T r ^ 103MeV [20 



By now imposing the condition M cr /2M 7T < 1 we have a constraint on the temperature 
ratio T/T c 



2ikL 



< 1 



T 

— > 0.95 

T r ~ 



(30) 
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and this constraint fixes T* = 0.95 T c ~ 98MeV. 

The approximations we have made up to now (we have evaluated (T) and b\(T) at T c , 
and we have chosen for sq(T) the behaviour as given by Eq. (|28"D) allow us to study M a jM. K 
only for T very close to T c . To study the mass ratio Mu/M^ as a function of T for all the 
values of the temperature from T = up to T very close to T c we have to go back to Eq. fl23"D 
which is valid for any T at the leading order in the quark mass m. To evaluate M%/M% we 
can perform the integration in yo by using the Cauchy theorem and summing up the series; 
then, through a numerical calculation we can finally perform the integration left. In Fig. [I] 
we plot the mass ratio M a jM. K as a function of T for a mass m = (m u + md)/2 = 9.5 MeV; 
from the numerical analysis we can evaluate again the temperature T* 

M T 

' < 1 => ^ > 0.93 (31) 



2M n ~ T c 

and this constraint now fixes T* = 0.93 T c ~ 95 MeF. 

We remark that this result for T* is in complete agreement with that previously derived 
by studying the Landau expansion of the effective potential for T very close to T c . Finally, 
if one wants to study M a /M n including all the orders in the quark mass m and for any 
temperature it is necessary to go back to the general formula fl22p . In Fig. ^] it is plotted 
the behaviour of the mass ratio as a function of T as derived by the explicit evaluation of 
Eq. (p2|) . In this case the temperature T* we are looking for is determined again by the 
condition 

M„ T 

< 1 => — > 0.95 (32) 



2M n ~ T c 

which fixes T* = 0.95 T c ~ 98 MeV. Curiously this result exactly coincides with that 
obtained through the Landau expansion and the evaluation of the coefficients a± (T), b\(T) 
at T c and by putting for s the expression (p8|) as suggested by the mean-field theory. 



13 



V. CONCLUSIONS 



We have applied the composite operator formalism to QCD at finite temperature. After 
carrying out consistently the renormalization in the effective potential, we have calculated 
the masses of the scalar and pseudoscalar mesons. The relevant QCD parameters have 
been derived from the expressions for the masses and decay constants of the pseudoscalar 
multiplet, by inserting as inputs the mass and decay constant of the pion, the charged kaon 
mass and the electromagnetic mass difference between neutral and charged kaon. 

A main purpose of the calculation was to determine the temperature T* at which 
M a (T) < 2M 7r (T). Such temperature will be reached before the critical point, since the 
scalar meson mass M a (T) is expected to decrease along with partial chiral symmetry restora- 
tion in hot matter, and at the same time the pion will lose its Goldstone nature. The temper- 
ature T* is also expected not to lie far from the critical temperature of the chiral transition. 
In principle the a-meson would appear as a narrow state at high enough temperature. For 
T > T* its partial width by strong interaction has to vanish. 

The result of our calculation gives T* = 0.95 T c , in an approximation including current 
quark masses and with the full effective potential given by the formalism. We find that this 
result almost coincides with a much simpler treatment using a mean-field expansion of the 
effective potential a-la Landau around the critical point in the limit of small quark mass. 

It is perhaps premature to discuss whether a possible experimental signature for the 
QCD phase transition may be inferred from the behaviour of the sigma width. The mass 
ratio itself is in any case a significant parameter to describe the degree of chiral symmetry 
breaking at different temperatures and therefore it has to be considered as a quantity of 
intrinsic theoretical interest. 
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FIGURES 

FIG. 1. Plot of M a /M 1T versus T at the leading order in m, for m = 9.5 MeV . 



FIG. 2. Plot of M a /Mir versus T at all orders in m, for m = 9.5 MeV. 
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